Fast rotation of an ultra-cold Bose gas 
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We study the rotation of a ^^Rb Bose-Einstein condensate confined in a magnetic trap to which a 
focused, off resonant gaussian laser beam is superimposed. The confining potential is well approxi- 
mated by the sum of a quadratic and a quartic term which allows to increase the rotation frequency 
of the gas above the trap frequency. In this fast rotation regime we observe a dramatic change in 
the appearance of the quantum gas. The vortices which were easily detectable for a slower rotation 
become much less visible, and their surface density is well below the value expected for this rotation 
frequency domain. We discuss some possible tracks to account for this effect. 
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The fast rotation of a macroscopic quantum object of- 
ten involves a dramatic change of its properties. For 
rotating liquid ^He, superfluidity is expected to disap- 
pear for large rotation frequencies A similar effect 
occurs for type-II superconductors placed in a magnetic 
field, which loose their superconductivity when the field 
exceeds a critical value For the last few years it 
has been possible to set gaseous Bose-Einstein conden- 
sates into rotation, by nucleating in the gas quantized 
vortices which arrange themselves in a triangular lattice 
H^QIEr^, 7]. The theoretical investigation of the prop- 
erties of these fast rotating gases has led to several possi- 
ble scenarios. Depending on the confinement of the gas, 
its dimensionality and the strength of the interactions, 
theoretical predictions involve the nucleation of multiply 
charged vortices |^ 0, 0, ^ melting of the vortex lat- 
tice 12, 13], or effects closely connected to Quantum Hall 
physics Q,EEIi3. 

The confinement of cold atomic gases is generally pro- 
vided by a magnetic trap which creates an isotropic trans- 
verse potential muj'j_r'^/2 in the plane perpendicular to 



the rotation axis z (we set + y 



2 _ ^2 



). As the rotation 



frequency Q varies from to lj± , the number of quantized 
vortices increases J^]. The limiting case of Q = u;±^ 
e.g. fast rotation, is singular: the confinement vanishes 
since the trapping force is compensated by the centrifu- 
gal force. In classical terms, the only remaining force on 
the particle is the Coriolis force. In quantum terms this 
corresponds to a gauge field A{r) = mft x r/2. This is 
analog to the physics of charged particles in a uniform 
magnetic field and one expects to recover for weak inter- 
actions an energy spectrum with Landau levels separated 
by 2 hij± . However the absence of confinement in a pure 
harmonic trap rotating at = u;± makes this study ex- 
perimentally delicate 19]. The Boulder group recently 
reached Q = 0.99 u;± using evaporative spin-up [2^. 

In the present paper, following suggestions in 
llH we study the rotation of a Bose-Einstein condensate 
in a trap whose potential is well approximated by a su- 
perposition of a quadratic and a small quartic potential: 
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We can thus explore with no restriction the domain of 
rotation frequencies Vt around uo^. Our results show a 
strong dependance of the behavior of the gas with Vt. For 
Q./uo±_ < 0.95 we observe a regular vortex lattice. This 
lattice gets disordered when ft increases. For Q > uj± the 
number of detectable vortices is dramatically reduced al- 
though we have clear evidences that the gas is still ultra- 
cold and in fast rotation. We conclude this Letter by 
proposing some possible explanations for this behavior. 

Our ^^Rb condensate is formed by radio-frequency 
evaporation in a combined magnetic+laser trap. The 
pure magnetic trap provides a harmonic confinement 
along the 3 directions with the frequencies a;^^/27r = 
75.5 Hz and uJz/^tt = 11.0 Hz. We superimpose a blue de- 
tuned laser (wavelength 532 nm) propagating along the z 
direction to provide the quartic term in the confinement. 
The potential created by the laser is 
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The laser's waist is = 25 /am and its power is 1.2 mW. 
The first term Uq (~ /c^ x 90 nK) in the right hand side of 
(O is a mere shift of the energy scale. The second term is 
a correction of the transverse trapping frequency; we get 
for the combined magnetic-laser trap u;±/27r = 65.6 Hz 
[2l|. The last term in Q corresponds to the desired 
quartic confinement, with k = 2.6 (3) x 10~^^ J.m~^. The 
expansion in is valid for r± < ^/2, which is indeed 
the case for all data presented here. 

We start the experimental sequence with a quasi-pure 
condensate that we stir using an additional laser beam, 
also propagating along z Q . This laser stirrer creates an 
anisotropic potential in the xy plane which rotates at a 
frequency l^stir- To bring the condensate in rotation at a 
frequency close to cjx, we use two stirring phases. First 
we choose — uj±/ a/2, so that the stirring laser reso- 
nantly excites the transverse quadrupole mode m = +2 
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FIG. 1: Pictures of the rotating gas taken along the rotation axis after 18 ms time of flight. We indicate in each picture the 
stirring frequency Q^^-^^ during the second stirring phase (a;x/27r = 65.6 Hz). The real vertical size of each image is 306 /xm. 



of the condensate at rest 22]. The duration of this first 
excitation is 300 ms and we then let the condensate relax 
for 400 ms in the axisymmetric trap. This procedure sets 
the condensate in rotation, with a vortex lattice contain- 
ing typically 15 vortices. 

We then apply the laser stirrer for a second time pe- 
riod, now at a rotation frequency close to or even 
above uo^. For the condensate with 15 vortices prepared 
during the first phase, the transverse quadrupole mode is 
shifted to a higher frequency and can now be resonantly 
excited by a stirrer rotating at ^^^2t ~ ^± • The 
second stirring phase lasts for 600 ms. It is followed by a 
500 ms period during which we let the condensate equi- 
librate again in the axisymmetric trap. At this stage the 
number of trapped atoms is 3 x 10^. 

We find that this double nucleation procedure is a re- 
liable way to produce large vortex arrays and to reach 
effective rotation frequencies l^eff around or above uo^. 
In particular it is much more efficient than the proce- 
dure which consists in applying permanently the stirring 
laser while increasing continuously its frequency. It has 
one drawback however: the effective rotation frequency 
r^eff of the gas after equilibration might differ significantly 
from the rotation frequency ^\^{^ that we apply during 
the second stirring phase. The role of the stirring phase 
is to inject angular momentum in the system, and the 
rotation frequency of the atom cloud may subsequently 
change if the atom distribution (hence the moment of 
inertia) evolves during the final equilibration phase. 

The rotating atom cloud is then probed destructively 
by switching off the confining potential, letting the cloud 
expand during r = 18 ms and performing absorption 
imaging. We take simultaneously two images of the atom 
cloud after time of flight (TOF). One imaging beam is 
parallel to the rotation axis z and the other one is per- 
pendicular to z [2^. 

Fig.^-h shows a series of images taken along the z axis 

f2) ('2) 

and obtained for various around uo^. When ^\^{^ in- 
creases, the increasing centrifugal potential weakens the 
transverse confinement. This leads to an increasing ra- 
dius in the xy plane. As it can be seen from the transverse 
and longitudinal density profiles of Fig. O obtained for 
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^stir/^'^ = 66 Hz, the gas after TOF has the shape of a 



in fast rotation is shown on the image (h) at the extreme 
right of Fig. n It corresponds to r2g^^^/27r = 69 Hz, far 
from the quadrupole resonance of the condensate after 
the first nucleation. In this case we could not bring the 
gas in the desired fast rotation regime. 

f2) 

For the pictures (e,f) obtained with l^^^-y 27r = 66 and 
67 Hz, the optical thickness of the cloud has a local 
minimum in the center; this indicates that the confin- 
ing potential in the rotating frame, Kot(^) = '^(^) ~ 
^^eff^l/^' has a Mexican hat shape. The effective rota- 
tion frequency l^eff thus exceeds uo^. The striking feature 
of these images is the small number of visible vortices 
which seems to conflict with a large value of (^eff- The 
main goal of the remaining part of this Letter is to pro- 
vide more information on this puzzling regime. 

In order to analyze quantitatively the pictures of Fig. ^ 
we need to model the evolution of the cloud during TOF. 
For a pure harmonic potential (/c = 0), a generalization 
of the analysis of [2^ to the case of a rotating conden- 
sate shows that the expansion in the xy plane is well 
described by a scalin g of the initial distribution by the 
factor (1 +cj^r^)^/^ [2^. We assume here that this is 
still approximately the case for a condensate prepared in 
a trap with a non zero quartic term kr\/A. 

We have analyzed a series of 60 images such as those 
of Fig. n assuming an initial atomic distribution given 
by the Thomas-Fermi law: 
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where a = 5.2 nm is the scattering length characteriz- 
ing the atomic interactions and ja the chemical potential. 
The optical thickness for the imaging beam propagating 
along proportional to the column atomic density in the 
xy plane, varies as (a + /3r^ + 7^1)^^^ 7 where a, /3, 7 can 
be expressed in terms of the physical parameters of the 
problem [23. An example of the fit, which takes /i and 
l^eff as adjustable parameters, is given in Fig. The 
agreement is correct, though not as good as for conden- 
sates confined in purely harmonic traps. This may be a 
consequence of the approximate character of the scaling 

transform that we use to describe the TOF evolution. 
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The resulting values for ^leff as a function of l^g^-^ are 
given in Fig. IHl We find l^eff — ^sti stirring frequen- 
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fiat pancake. The limit of our method for setting the gas cies < 68 Hz. Note that for ^l^U^^ = (image 
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FIG. 2: Optical thickness of the atom cloud after time-of- 
flight for r2g^-^^/27r = 66 Hz. (a) Radial distribution in the xy 
plane of Fig.^. Continuous line: fit using the Thomas- Fermi 
distribution (b) Distribution along the z axis averaged 
over \x\ < 20 /im (imaging beam propagating along y). 



Fig. ^) the quality of the fit is comparatively poor due 
to local inhomogeneities of the atom cloud. 

From the value of /i given by the fit, we recover the 
atom number J^J. For the largest measured rotation fre- 
quency, the Thomas-Fermi distribution Q corresponds 
to a nearly spherical atom cloud before TOF (diameter 
29 jam in the xy plane and length 34 jam along z). 

A second determination of the effective rotation fre- 
quency l^eff of the condensate is provided by the vor- 
tex surface density after TOF. Assuming that the vortex 
pattern is scaled by the same factor as the condensate 
density, we deduce the vortex density py before TOF, 
hence the rotation frequency l^eff = Trhpy/m 18]. For 
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^stir ^ value of l^eff dcduccd in this way and 



plotted in Fig. |31 is in fair agreement with the one de- 
duced from the fit of the images. On the contrary, for 
r^g^-^/27r = 66 — 68 Hz the number of distinguishable vor- 
tices is much too low to account for the rotation frequen- 
cies determined from the fits of the TOF images. 

In order to gather more information on the rotational 
properties of the gas, we now study the two transverse 
quadrupole modes m = ±2 of the gas. We recall that 
these two modes have the same frequency for a non ro- 
tating gas, due to symmetry. For a rotating gas in the 
hydrodynamic regime, the frequency difference — 
is proportional to the ratio between the average angular 
momentum of the gas and its moment of inertia, which is 
nothing but the desired effective rotation frequency (^eff 
0. The strength of this approach is that it does not 
make any assumption on the expansion during TOF. 

To study these transverse quadrupole modes, we briefly 
illuminate the rotating gas using the laser stirrer, now 
with fixed axes [2^. This excites a superposition of the 
modes m = ±2 with equal amplitudes. We then let the 
cloud evolve freely in the trap for an adjustable duration 
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FIG. 3: Effective rotation frequency Qeff as a function of the 
stirring frequency ^^^til- ^- values deduced from the fit using 
the Thomas- Fermi distribution O: values obtained by 

measuring the vortex density in the TOF pictures. A: values 
obtained using surface wave spectroscopy. The bars indicate 
standard deviations. 



and we perform the TOF analysis. From the time varia- 
tion of both the ellipticity of the cloud in the xy plane and 
the inclination of its eigenaxes, we deduce cj±2, hence the 
effective rotation frequency l^eff = (^+2 — ^-2)/2. The 
corresponding results are plotted in Fig. IHl They are in 
good agreement with the results obtained from the fits of 
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the images: l^eff — ^st/r- Consequently they conflict with 
the rotation frequency that one derives from the vortex 

(2) 

surface density when r^g^/./27r = 66 — 68 Hz. 

To explain the absence of visible vortices in the regime 
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^stir/^'^ = 66 — 68 Hz, one could argue that the gas might 
be relatively hot, hence described by classical physics. 
The equilibrium state should then correspond to rigid 
body rotation. However all experiments shown here are 
performed in presence of radio- frequency (rf) evapora- 
tive cooling. The rf is set 24 kHz above the value which 
empties the trap and it eliminates all atoms crossing the 
horizontal plane located at a distance Xev ~ 19 pm be- 
low the center (one-dimensional evaporation). /,From 
evaporative cooling theory 29, 30] we know that the 
equilibrium temperature T of the rotating gas is a frac- 
tion of the evaporation threshold Kot(^ev) 32 nK for 
l]eff/27r = 67 Hz 31]. This indicates that T is in the 
range of 5 — 15 nK, well below the critical tempera- 
ture at this density (180 nK for an estimated density 
3 X 10^^ cm~^) |33|. This low temperature is also con- 
firmed by the small decay rates of the quadrupole modes 
(- 20 s-^), characteristic of T <C T^. 

The estimated temperature T is of the order of the 
splitting between Landau levels 2huj_\_/kB = 6.3 nK, so 
that only the first two or three levels are appreciably 
populated. For each Landau level, the rf evaporation 
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eliminates states with an angular momentum Lz/h > 
(xev/<^ho)^ ^ 200, where aho = \/^J{Tfujo^ is the ground 
state's size of the harmonic oscillator with frequency uj^. 

A first possibility to interpret our data consists in as- 
suming that the whole gas shown in Fig.QJ-g is rotating, 
that it is in the degenerate regime, and that it cannot 
be described by a single macroscopic wave function since 
the number of distinguishable vortices is too low to ac- 
count for the measured rotation frequency. In this point 
of view, the state of the system is therefore more complex 
than a single Hartree state. A theoretical analysis along 
this line, involving the formation of composite bosons, 
has been very recently proposed in js^ j. 

A second possible explanation is that the vortices are 
still present for > cj^, but that the vortex lines are 
strongly tilted or bent, so that they do not appear as 
clear density dips in the images of Fig. ^-g. This could 
be for example a consequence of the residual anisotropy 
of the trapping potential in the xy plane 21]. In an ex- 
periment performed at Q = 0.95 cj^, the Boulder group 
has shown that a static anisotropy of ~ 4% reduces con- 
siderably the vortex visibility over a time scale of 0.8 s 
|3^ . However the residual anisotropy of our trap is < 1% 
and this effect should be limited, unless the vortex lat- 
tice becomes notably more fragile as ^} approaches u±. 
The distortion of the vortex lattice could also result from 
a strong increase of the cristallization time for Q ~ lj± 
which would then exceed the 500 ms equilibration time 
Teq. To check this effect, we have increased Teq up to 2 s 
without noticing any qualitative change in the atom dis- 
tribution (for larger Teq the gas slows down significantly). 

A third possibility is that a small fraction of the gas has 
stopped rotating during the standard 500 ms equilibra- 
tion time, and can therefore be at a higher temperature 
1^3 • This small fraction could decrease the visibility of 
the vortices present in the condensed part of the cloud. 
The existence of such a fraction could be tested by a 
quantitative analysis of the density profiles of Fig. [S^b. 
Such an analysis should take into account the TOF evo- 
lution of an interacting atomic cloud initially confined in 
a non harmonic potential. 

To summarize we have presented in this Letter a direct 
evidence for a qualitative change in the nature of a degen- 
erate Bose gas when it is rotated around and above the 
trapping frequency. We plan to complement this study 
by an analysis of other modes of the rotating gas, such 
as the transverse breathing mode |35(. Several other in- 
teresting phenomena have been predicted and remain to 
be investigated experimentally for this type of potential: 
existence of several phases involving either a vortex array 
with a hole, or a giant vortex |3. I^ITol. Illl . Another 
extension of the present work consists in transposing the 
experimental scheme to a 2D geometry, where the motion 
along z would be frozen. The situation would then be the 
bosonic analog of the situation leading to the quantum 
Hall effect, provided the effective rotation frequency of 



the gas is precisely adjusted to the trap frequency. 
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